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Abstract 

Context. Many protoplanetary nebulae show strong asymmetries in their surrounding shell, pointing to asymmetries during the mass 
loss phase. Questions concerning the origin and the onset of deviations from spherical symmetry are important for our understanding 
of the evolution of these objects. Here we focus on the circumstellar shell of the post-AGB star HD 161796. 
Aims. We aim at detecting signatures of an aspherical outflow, as well as to derive the properties of it. 

Methods. We use the imaging polarimeter ExPo (the extreme polarimeter), a visitor instrument at the William Herschel Telescope, to 
accurately image the dust shell surrounding HD 161796 in various wavelength filters. Imaging polarimetry allows us to separate the 
faint, polarized, light from circumstellar material from the bright, unpolarized, light from the central star. 

Results. The shell around HD 161796 is highly aspherical. A clear signature of an equatorial density enhancement can be seen. 
This structure is optically thick at short wavelengths and changes its appearance to optically thin at longer wavelengths. In the 
classification of the two diff'erent appearances of planetary nebulae from HST images it changes from being classified as DUPLEX at 
short wavelengths to SOLE at longer wavelengths. This strengthens the interpretation that these two appearances are manifestations 
of the same physical structure. Furthermore, we find that the central star is hotter than often assumed and the relatively high observed 
reddening is due to circumstellar rather than interstellar extinction. 
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1. Introduction 

When image d, post-AGB stars often display a highly nonspheri- 
cal structure ( |Ueta et aL||2000| ). It has been, and still is, a matter 
of debate what is the cause for this asphericity and at which point 
in the post main sequence evolution the deviations from a spher- 
ical wind start. One of the most favored explanations is through 
interaction of the wind with a binary companion. However, re- 
cently Hriv nak et al. ( 2011 1 have searched for such compan- 
ions and found only one out of seven proto-planetary nebula 
contained a detectable binary companion. Detailed images of 
evolved stars together with extended modeling is required to get 
a better handle on the origin of the non- spherical nature of many 
(proto)-planetary nabula. 

The types of assymetries seen in protoplanetary nebulae 
have been devided into two classes by |Ueta et aL] (|2000 ); the 
SOLE (Star-Obvious Low-level-Elongated) and the DUPLEX 
(DUst-Prominent Longitudinally-EXtended) types. The SOLE 
type shows one elongated structure which was interpreted as 
the appearance of an optically thin torus surrounding the star. 
The DUPLEX type on the other hand shows two distinct blobs 
at opposite sides of the central star. This was interpreted as 
the appearance of a bipolar outflow structure. In the DUPLEX 
type there is also a torus, but it is optically thick and thus self- 
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shielding and not visible in scattered light. Essentially, both 
types are appearances of the same physical structure, where the 
main difl'erence lies in the optical depth through the torus. 

Here we present polarimetric images of the post-AGB object 
HD 161796 observed with the extreme polarimeter (ExPo) at the 
William Herschel Telescop e (WHT). HD 161796 is an F3Ib star 
( |Fernie & Garrison] |1984| ) with a detached circumstellar shell 
as deri ved from its double peaked spectral energy distribution 
(SED) ( [Hrivnak et al. 1989| ). The circumstellar matter is clearly 
observed to be oxygen rich and dominated by silicate emission 
(Justtanont et al.| |1992t [Hoogzaad et al. l 120021). Spat ially re- 
solved images in the ne ar-IR (|Ueta et air , 2000 ; Gled hill et al. 



2001 



, and the mid-IR fSkinner et aLnT994, .Gledhill & Yates 



2003 ) resolve the inner boundary of the shell to be much larger 
than the dust condensation radius. From these observations it can 
be derived that the mass loss stopped around 430 years ago. 

The shell surrounding HD 161796 has been modeled from 
difl'erent perspectives. Hoogza ad et al.| ([2002) provide an ex- 
tensive analysis of the mineralogy and size distribution of the 
dust grains in the envelope using a spherical symmetric radiative 
transfer model. They derive from their fit that the mass loss rate, 
which ended 430 years ago, was realtively high (5 • 10"^MQ/yr) 
which they suggest is consistent with the prominent presence of 
ice-coated grains. Deviating from spherical symmetry, Meixner] 
et al.| ( |2002] ) developed a general, axisymmetric parameterization 
to fit images of protoplanetary nebulae. They successfully apply 
this to HD 161796 and manage to simultaneously fit the SED 



and the mid IR images. The axisymmetric model by [Meixner 
eTaL] ( [20021 ) is very successful in reproducing the various mor- 
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phologies observed. However, as was shown by Gledhill & Yates 
2QQ3| ), there are also indications of deviation from axisymmetry 



Table 1. Details of the observations 



because in the mid-IR images the two lobes east and west from 
the central star are not equally bright. 

We observed HD 161796 using imaging polarimetry with 
various filters at optical wavelengths. The observed morphology 
of the circumstellar shell is observed to change with wavelength 
(section [2]). In this paper we present a model which is consistent 
with the images at diff'erent wavelengths and fits the SED from 
optical to millimeter wavelengths (section [5]). The wavelength 
dependent morphology can be understood as a changing optical 
depth through the equator of the system. This is used to accu- 
rately derive the total optical depth through the torus and con- 
strain the scattering properties of the circumstellar dust grains. 



2. Polarimetric imaging using ExPo 

Imaging polarimetry takes advantage of the fact that scattered 
light from circumstellar material is polarized to filter out the 
bright point spread function (PSF) of the central star. This way 
we can image the surroundin g material down to a contrast of ~ 
10""^ ( jRodenhuis et al. 2008| ). For our observations, the imaging 
polarimeter ExPo was mounted at the 4.2 m William Herschel 
Telescope (WHT). The raw frames from the observations were 
carefully reduced using a state-of-the-art data reduction pipeline 
including detailed recentering of the frames to reduce the eff'ects 
of seeing. The technical setup and data reduction steps are briefly 
outlined below. For details we refer to |Canovas et al.| ( |2011j . 



Date 



Filter 



2.1. Technical description 

ExPo is a regular visitor instrument at the William Herschel 
Telescope, mounted on the Nasmyth platform (Rodenhuis et al. 



|2008 |). ExPo is a dual-beam imaging polarimeter, operating at 
visible wavelengths 500 - 900 nm, that combines fast modula- 
tion with the dual-beam technique to minimize systematic errors. 
The design concept of ExPo is that a Ferroelectric Liquid Crystal 
(FLC) modulates the light, rotating its polarization plane by 90 
degrees. The FLC changes its state every 0.028 seconds (approx 
35 HZ), which allows us to reduce the atmospheric (seeing) ef- 
fects. A cube beamsplitter divides the light into two beams with 
orthogonal polarization states. These two simultaneous images 
are then projected onto an Electron Multiplying Charge Coupled 
(EM-CCD) camera. 

2.2. Observations 

The observations of HD 161796 were taken in May 2010. To 
determine the wavelength dependence of the morphology of 
HD 161796 observations were also secured using a set of dif- 
ferent filters. For both the H^, and the Na filters we took images 
using a filter at the spectral position of the line and in the contin- 
uum next to it. No significant diff'erences were detected between 
the line and continuum observations, so we combined them to 
increase the signal to noise ratio. Details of the observations and 
filters used are shown in Table [T] To obtain a linear polarization 
image in each filter we take four sub-exposures with FLC angles 
of 0°, 22.5°, 45° and 67.5°. For HD 161796, each FLC angle 
comprises at least 4000 frames of 0.028s duration. Additionally, 
a set of calibration flat-fields is taken at the beginning and end of 
every night and a set of dark frames is taken at the end of each 
observation. 



25 May 2010 Na continuum (0.590+0.001 jim) 

25 May 2010 Na line (0.580+0.001 fim) 

24 May 2010 continuum (0.670+0.001 yt/m) 

24 May 2010 line (0.660+0.001 fim) 

25 May 2010 Sloan I (0.76+0.08 //m) 



2.3. Data analysis 

The FLC modulates the polarization state of the incoming light 
by 90° every 0.028s switching between "A" and "B" frames. For 
each frame ("A" or "B") the light is separated, using a beamsplit- 
ter, into two orthogonal beams which are simultaneously imaged 
on the EM-CCD, a concept that is the basis of a dual-beam po- 
larimeter. For every pair of observations this results in four dif- 
ferent images (i.e. in total at least 8000 for each FLC angle as 
described above). This results in four diff'erent images per FLC 
cycle ("A"-h"B"): ^left^^right^ ^left^ ^right- Each frame is first flat- 
field corrected, dark subtracted, and cleaned of cosmic rays. The 
four images are combined using the double-diff'erence approach 
dKuhn et al.|[2QQT||Hinkley et al.|[2QQ9l ). 



= 0.5(AA - AB) = 0.5((Aieft - A^ght) - (^Mt - bright)) (1) 

where P'^ is the uncalibrated linearly polarised image. We use 
here P'j in stead of Q or U since the rotation of the sky and the 
imperfect alignment of the instrument cause the images obtained 
in this way to be neither exactly Stokes Q or Stokes U. These 
images are calibrated to produce the Stokes Q and U images 
in the reference frame attached to the target. The total intensity 
image is computed as: 



/ = 0.5(Aieft + Aright + ^left + ^rig] 



iht) 



(2) 



The image alignment algorithm is based on a cross-correlation 
algorithm that aligns each image with a template, with an accu- 
racy of a third of a pixel. The sky polarization is subtracted from 
the images by computing the polarization of four difl'erent sky 
regions o n the images. The d ata analysis is described in more 
detail by Canovas et al.| ( 20lT] ). 

Finally, the reduced images are calibrated using the method 
of Rodenhuis et al. (in prep) in order to produce Stokes Q and U 
images. The polarised intensity is defined 
and the degree of polarisation as P = Pi/ 1. The polarisation 
angle, which defines the orientation of the polarisation plane is 
Pq = 0.5 X arctan U/Q. The instrumental design of ExPo in- 
cludes a polarisation compensator which reduces the instrumen- 
tal polarization from ~3% to approximately 1%. This is further 
removed in the data analysis by assuming that the central star is 
unpolarized. We stress here that the fact that we remove the po- 
larization of the central resolution element causes any polariza- 
tion present inside the inner resolution element to be removed. 
This also implies that the polarized intensity we derive over the 
entire image has to be take n as a lower limit. This eff'ect is de- 
scribed in Mi n et al. 1(20121 ). 



2.4. Polarization images 

The final polarized intensity images obtained with ExPo for 
three diff'erent filters are displayed in the upper panels of Fig. T] 
The wavelength increases from left to right. In these images the 
inner edge of the shell is clearly resolved. The inner edge of the 
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Figure 1. Polarized intensity images of HD 16 1796 in the different filters (upper panels) as obtained with ExPo. Left panel shows 
the ExPo observations using the Na 580 nm filter, middle panel the Hq, 656 nm filter, and the right panels the Sloan I filter. The 
lower panels display the best fit model images. The absolute scaling of the model images is adjusted to provide a good fit to the 
observations. 



shell is derived from our images to be at ~ 0.8'', which translates 3. Modeling the circumstellar dUSt shell 



to 960 AU at a distance of 1200 pc ( Skinner et al.|t 199 4 ). We find 
a bright ring with an intensity drop which is quite sharp outside 
of 1.5'' (1800 AU). Note that the image taken in the Na filter has 
a seeing that is significantly worse than the other observations. 
At the shortest wavelength, the Na filter at 580 nm, the image 
displays a brightness enhancement towards the north and south. 
This structure is consistently seen in both the Na line and the 
Na continuum filters. The brightness enhancement is not clearly 
visible in the other two images. The 656 nm image shows still 
some structure which is consistently seen in the line and contin- 
uum image. The Southern part of this image seems to be system- 
atically slightly weaker than the Northern part. Also, the Sloan I 
image still displays some very low level fluctuations in the ring, 
with a brighter arc to the top and a small brigthness enhancement 
to the lower right part of the image. Furthermore, the structure 
seen in all images is elongated in north-south direction with an 
aspect ratio of ~ 1.3. 



The procedure to correct for instrumental polarization alters 
the inner region of the images. In order to verify that the ob- 
served ring structure with an empty inner region is not an efl'ect 
of this procedure, we have also looked at the images before in- 
strumental polarization correction. We consistently see the ring 
structure in all images also before the correction. Furthermore, 
the derived inner edge of the ring, as well as the structure seen 
in it, is unafl'ected by the removal of instrumental polarization. 



In this section we first describe the parameterization chosen for 
the density setup and the assumptions for the optical properties 
of the dust grains. Second, we also discuss the radiative transfer 
model we use. 



3A. Density setup 



We take the approach from Meixner et al. (2002) who find a 
model which is consistent with the thermal infrared emission 
(both imaging and the spectral energy distribution). This is a pa- 
rameterized description of a bipolar structure with an equatorial 
density enhancement and is given by 



Po 



B< 1+Csin^ 



R 



x^l -hA(l -cos^)^ 



-{RminlRs^Y 



(3) 



were R is the radius to the star, Rs^ is the radius of the super- 
wind phase, Rmin is the inner radius of the dust shell, 6 is the 
angle measured from the pole, po is the density at the inner edge 
of the shell, and A-F are fitting parameters. This density setup is 
empirical, and thus the parameters do not directly relate to phys- 
ical quantities. The density model from Meixner et al.| (2002) 
assumes an equatorial density enhancement caused by the super- 
wind in a more or less torus like structure. In the remaining of 
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this paper we will use the term 'torus' to refer to this equatorial 
density enhancement. Outside of the superwind phase these au- 
thors place a constant mass loss. As discussed below, we adopt 
their model, and only tune it slightly to account for the differ- 
ences in dust grain properties we use. A cartoon drawing of the 
density structure is shown in Fig. [2] 



large fraction of this reddening is circumstellar rather than inter- 
stellar. Therefore, we will assume the usual distance of 1.2 kpc. 
The luminosity of the star was varied in order to fit the SED. By 
varying the eff'ective temperature of the star we find that indeed 
the best fit is obtained by using a Kuruzc model with eff'ective 
temperature 7250 K, logg = 0.5, and a luminosity of 3000 Lq. 



3.2. Dust grain properties 



3.4. Radiative transfer 



We have computed the absorption and scattering properties of 
grains with a mixture of silicates, metallic iron and water ice 
with mass abundances of 77, 15, and 8 % respectively. We do 
not aim here for an exact compositional analysis, for that we re- 
fer to the thorough analysis of the ISO spectrum by 'Hoogzaad] 



For the radi ative transfer we use the highly flexible MCMax code 
( |Min et alTl|2009) . Thi s code is based on the Mon te Carlo radia- 
tive transfer method by |Bjorkman & Wood|p001| ) and takes into 
account full anisotropic scattering and polarization. The code 



et al.| ([2002 ). For the composition of the silicates we take the 



mixture derived by , Min et aL] ( |2007p for interste llar silicates. We 
use the refractive indices for the materials from lDorschner et aTl 
1995|) a nd'Henn ing & Sto gnienko'('1996') for the silicates, from 



IQrdal et al. (1988') for the metallic iron, and from Warren (1984 ) 
for the water ice. The exact grain composition has little influence 
on the scattering properties of the particles compared to the par- 
ticle size and shape. However, the composition does have a large 
impact on the spectral energy distribution (SED). In general, the 
size and shape distribution of the grains is well constrained by 
the polarimetric measurements, while the composition is much 
better constrained from the spectral features seen in the infrared. 

In many studies it is assumed that the grains are homoge- 
neous spheres. However, although applicable in some particu- 
lar cases, the homogeneous dust grain model generally produces 
resonances that are not observed in the optical properties of nat- 
ural particles (Min et ar,"2012 ). Though one can often remove 
the sharp nature of these resonances by using a size distribu- 
tion of spheres, the smoothed resonances still cause efl'ects that 
are inconsistent with computations or measurements for irreg- 
ularly shaped particles, this can lead to errors in the interpre- 
tation of observati ons where t he dust grains are expected to be 
irregularly shaped. |Min et alT] (|2005 ) have shown that breaking 
the perfect shape of a homogeneous sphere even in the simplest 
way already reduces the impact of these resonances significantly. 
Therefore, for this study we will use the DHS grain shape model, 
as suggested by |Min et al.| ( |2QQ5| ), using an irregularity parame- 
ter/max = 0.8. 

We take a grain size distribution with minimum grain size 
0.005 yum and maximum grain size 0.5 yum. We take a powerlaw 
size distribution and vary the index to match the images. 



3.3. Central star 

The eff'ective temperature and luminosity of the central star 
are a matter of debate and varies from 6666 K ( Cenarro et al.[ 
2007| ), derived from optical spectroscopy of the stellar spectrum, 
to 7500 K (Gle dhiU & Yatesl [2003|, derived from fitting the 
SED. The mos t recent detailed analysis of the optical spectrum 
( Kipper| 2007) results in an effective temperature of 7250 K, 
which is at the high end of the temperatures found in the lit- 
erature. In the same paper these authors also claim that the lumi- 
nosity of the star is significantly higher than previously assumed, 
and from this the derived distance must be larger as well. In or- 
der to arrive at this claim |Kipper| ( [2007] ) assume a relatively high 
initial mass of the central star (0.9 Mq ), which is at th e upper 
most end of the mass range derived by |Stasinska et al.| (2006). 
Furthermore, |Kipper| ( 2007 ) argue that the large reddening of the 
star points to a large distance, while we will show below that a 



has been ex tensively benchmarked against other codes ( [Pinte 
et al.||200^ . The polarimetric images are produced by integrat- 
ing the formal solution of radiative transfer after the local ra- 
diation field at each lo cation in the dust shell is found. Further 
details can be found in |Min et al.| ( |2QQ9l|2QT2l ). 

3.5. Instrument simulation 

An important part of comparing the resulting model images to 
the observed images is a proper simulation of instrumental and 
data reduction effects. We use a detailed simulation of seeing, 
noise and telescope effects to reduce our high resolution images 
to the spatial scale and noise level of the observations. The model 
images are convolved with a PSF matching the seeing conditions 
at the time of observing. Also, photon noise was added match- 
ing the computed fluxes. After this we apply the data reduction 
pipeline to the simulated observations in exactly the same way as 
was done for the ExPo observations. The instrument simulator is 



described in more detail by |Min et al.| ( |2012l ). 



4. Resulting best fit model 

Since we have a large, multi wavelength observational data set, 
and all observations are of different types, it is very hard to define 
the best fit in a proper, statistical sense. For this reason we have 
chosen to start off by using a by-eye fitting procedure, where we 
vary the parameters of the model until we find a solution that 
best reproduces the observations. This method does not give ro- 
bust estimates of the uncertainties of the derived parameters, and 
might be sensitive to local minima in parameter space. The ad- 
vantage however is that we can judge the goodness of fit for cer- 
tain datasets and in this way provide relatively robust constraints 
in certain wavelength intervals. For example, we conclude that 
we cannot reproduce the mid infrared images properly. We will 
come back to this issue later in the paper and discuss the impli- 
cations for the model. 

We start from the model by IMeixner et al.| (|2002j) and var- 
ied the parameters for the dust grain size distribution from there 
until we get a correct match for the polarimetric images. We 
also checked if we could get a better match by changing the 
density structure parameters, but find that, by eye, a good fit 
was obtained by taking the parameters A, B, C, E, and F from 



Mei xner et al.| ( |2002| . We take the parameter D to be equal to E, 
which gives a slightly different density structure than obtained 
by [Meixner et al.| ( 2002| ). If we comput e the mass loss history us- 
ing the value of D given by 'Meixner et al.| ( |200 2), which is D=l, 
we get a strange dip in the mass loss rate right before the high 
mass loss phase sets in. Putting D=E, like we did, avoids this 
problem (see also section [5]). Since the region affected by this is 
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right behind the optically thick torus, it makes almost negligible 
difference for fitting the SED and the ExPo images. 

In order to more robustly fine-tune the parameters of the 
model, the second step is to use a genetic fitting algorithm to find 
the best fit model to the SED, ISO spectrum and ExPo images. A 
genetic fitting algorithm uses the concept of evolution to find an 
optimum solution. The procedure starts out with a population of 
random models. The parameters of the models in this population 
are treated as their 'genes'. The best models of this population 
are given the highest change to reproduce their genes (parameter 
values) into the next generation (survival of the fittest). We use 
a population of 24 models and evolve this through 100 genera- 
tions. By studying the resulting 2400 models, we can also get a 
better feeling which parameters are constrained by the observa- 
tions and which are more degenerate. For m ore details on genetic 
fitting algorithms see Charbonneaul (|1995|). In order to run this 
procedure one needs to define what is a good fit. With only a sin- 
gle observational data set this is simply the reduced However, 
since we deal with diff'erent datasets which all have diff'erent un- 
certainties and what one would call a good fit is diff'erent for dif- 
ferent datasets (images are usually harder to reproduce exactly 
than an SED). Therefore, we have to define the weight we give to 
each dataset (images and SED). We checked that taking diff'erent 
values for these weights does not make a significant diff'erent to 
the best fit obtained. We fit the radial and azimuthal profiles of 
the polarized intensity images. The best fit we find does deviate 
slightly from the fit we get by eye, but all important parameters 
are very similar. The parameters of this final best fit model are 
presented in the next section. 

The final model was computed such that it simultaneously 
reproduces the overall SED and the ExPo polarization images. 
Also we aim to be roughly cons istent with the near infrared 
imaging polarimetry conducted b y|Gledhill et al. l('2001) and the 
mid infrared intensity images by|Gledhill & Yates (2003 ). The 
resulting SED is shown in Fig. [3] and the simulated images from 
the model are shown in the lower panels of Fig. [T] Although 
we cannot do absolute photometry, we can compare the degree 
of polarization to some extend. We have to be careful since the 
computed degree of polarization is basically the polarized in- 
tensity, caused by scattering off' dust grains, divided by the total 
intensity, which is mostly the tail of the PSF from the central 
star. Also, the low flux regions in the ExPo images are domi- 
nated by noise and unpolarized background emission. Therefore, 
the computed degree of polarization is influenced by seeing and 
instrumental efl'ects and cannot be used directly. Since our in- 
strument simulator does take most of these efl'ects into account, 
we can do a rough comparison. By doing this we find that we 
systematically over predict the degree of polarization with our 
model by a factor 3 to 5. This is most likely due to the relatively 
high degree of polarization of our model particles. Particles with 
more natural shapes and surface structure usually show a degree 
of polarization in the range 10-20 %, while our model particles 
have ~ 50 %. We conclude that likely small scale irregularity of 
the dust grains is the cause for the low degree of polarization 
observed (see also Section |4~2| for more detailed discussion). 

A notable diff'erence between the model images and the ob- 
served images is the lack of elongation at longer wavelengths. 
We find that, for the model, when the torus becomes optically 
thin the image becomes roughly circular, while the observations 
still display an elongated structure. This could indicate that the 
torus actually starts a bit closer to the star than the bipolar out- 
ffow indicating a diff'erence in outffow velocity; the polar outffow 
being faster by a small fraction. We consider it beyond the scope 
of our parameterized model to explore this further. 



bipolar 
structure 



spherical 
outflow 
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Figure 2. Cartoon drawing of the density structure around 
HD 161796. Note that the appearance of the bipolar structure 
is caused by the fact that the photons escape through the polar 
regions, illuminating the low density material there. 
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Figure 3. The spectral energy dis tribution together w ith th e 
model fit. ISO spectrum fro m |Hoogzaad et al. (2002). 



Photometry taken from Fernie (1983); Hrivnak et al. 



Walmsley'etaL] ( [r991J ; ,Meixner et al., ( .1999^ ; |Ueta et al.| ( |2QQQ 



1989) 



In the model SED we see a rather strong 3 yum absorption 
feature caused by the water ice in the torus. This feature is not 



seen in the ISO spectrum. According to Hoogzaad et al. (2002) 
this is caused by an anomaly in the measured spectrum. 

4.1. Derived parameters 

The best fit model parameters are given in Table [2] We find 
that the best fit to the SED and observed images is obtained for 
a nearly edge on geometry, where we look through the torus. 
This is consistent with previous modeling of infrared images 



(Meixner etal. 2002 



rameters for the shel 



Gledhill & ^^^tes||2Q03| . The resulting pa- 



structure are slightly diff'erent than previ- 
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Table 2. Values for the best fit model parameters 



Parameter 




Symbol 


Value 


Fixed parameters 


Effective temperature 




7250 K 


Radius of the star 






35 Rq 


Luminosity of the star 




3000 Lo 


Distance to the star 




D 


1200pc 


Inclination angle 




i 


80° 


Minimum grain size 




^min 


0.005 idm 


Maximum grain size 




^max 


0.5 jim 


Irregularity parameter 


fmax 


0.8 


Fit parameters 


Inner radius 




R-min 


904 AU 


Superwind radius 




R-sw 


970 AU 


Outer radius 




Rout 


1 16500 AU 


Dust density at inner edge 


Po 


6.7 • 10-21 g/cm^ 


Size distribution powerlaw 


P 


3.7 


Density distribution 




A 


12.2 (8) 


parameters 




B 


2.1 (2) 


(in parenthesis are the values 


C 


2.9 (2) 


from Meixner et al. ( 


20021) 


D 


6.4 (1) 






E 


6.4 (3) 






F 


1.0 (1) 







\ wavelength range ^ 


wavelength range 






\ ^ 








1 1 1 1 1 1 1 1 1 1 





0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 

wavelength [|im] 

Fi gure 4. The total integrated degree of polarization as observed 
by Parthasarathy et al. (2005]) together with the computed values 
from the best fit model. The arrows in the cartoon drawings left 
and right of the dashed line give the orientation of the polariza- 
tion vector in the model computations. The vertical dashed line 
indicates the zero integrated polarization wavelength. 



ously reported. This is caused by the additional constraints from 
the polarimetric images. Although a reasonable fit can be ob- 



tained by using the density distribution parameters from Meixner 
et al.| ( [2QQ2| ), a better fit to the images is obtained by us i ng the 



parameters in Table [2] The values from [Meixner e t al.| ( |2002[ ) 
are indicated in brackets for the density distribution parameters. 
Especially the D and E parameters are much higher. This is to 
create the sharp outer edge of the ring detected in our polarized 
intensity images. The outer radius of the shell is not constrained 
by the observations as long as it is large enough to reproduce the 
far infrared emission. 

The powerlaw of the size distribution is well constrained. We 
have to vary the grain properties in order to have significant opti- 
cal depth through the equator. Since the absorption optical depth 
is well constrained by the total amount of reprocessed energy in 
the infrared, we have to vary the albedo to create additional ex- 
tinction through scattering. We find that we need grains with a 
relatively high albedo to get sufficient optical depth at 0.58 jim 
to reproduce the observed morphology. 

The model we obtain is not unique. Many parameters in our 
model are to some degree degenerate. However, there are certain 
things we can constrain with high certainty. 

Optical depth through the equator: From the ExPo images 
we can constrain the optical depth through the equator quite ac- 
curately. We find that the models can only reproduce the equa- 
torial brightness gap seen in the Na image when the total (i.e. 
scattering -h absorption), optical depth is r ~ 3. 

Albedo of the grains: Since we know the total energy repro- 
cessed by the dust grains in the infrared, we actually have good 
constraints on the optical depth for absorption as well. The re- 
processed energy in the infrared should be equal to the amount 
of energy absorbed in the visual. Together with the total opti- 
cal depth, this gives us the single scattering albedo of the grains, 
which in our model co = Csca/Qxt ~ 0-84 at A = 0.58yum. 
This is actually quite a high value for small cosmic dust grains. 
Given the assumptions of the composition and shape distribution 
of the particles that we take, the albedo of the grains accurately 
constrains the size distribution. 



Inner radius of the outflow: We directly image the inner ra- 
dius of the torus. Also from previously published images, the 
inner radius of the torus is well constrained. We find it to be 
0.8'' which translates to 960 AU at the assumed distance. This 
inner radius scales linearly with the distance to the source. Note 
that the best fit model puts the inner radius slightly more inwards 
(904 AU). However this diff'erence is within the errors caused by 
our limited spatial resolution. 

4.2. Degree of polarization 



[Parthasarathy et aL] ( |2005| ) have measured the spatially inte- 
grated polarization of HD 161796 to be a bit less than 1% at 
optical wavelengths. The orientation of the polarization vectors 
was found to be roughly east- west (position angle of ~ 110°). 
We cannot compare these observations directly to the model cal- 
culations since the observations are likely heavily polluted by 
interstellar polarization. Since ExPo is mounted in Nasmyth fo- 
cus, we have to deal with a relatively high degree of instrumental 
polarization. This is removed by assuming that the central star is 
unpolarized, a method which is great for enhancing spatial vari- 
ations in the polarized intensity. However, it implies that we can 
not do absolute polarimetry using the ExPo images, and thus 
cannot compare our values with those from Partha sarathy et al] 
( i2QQ5| ). In Fig.|4]we show the values for the integrated polariza- 
tion predicted by our best fit model. We predict from the model 
that the angle of polarization will flip 90 degrees when going 
to the near infrared. This is because at the short wavelengths 
scattered light from the material ejected at the poles dominates, 
while at the longer wavelengths the scattered light from the torus 
takes over (see also the simulated near infrared images in Fig. [7]). 
The turnover point in our model is approximately when the opti- 
cal depth through the torus equals unity (see Fig. [5]). This result 
is consistent with the measure d polarization angle in the optical 
by Parthasarathy et al. ( 2Q05|) and with the po larization images 
in the near infrared from'G ledhill et al.| ( |2Q01| ). 

There are two solutions to the high degree of polarization 
we predict with respect to the observations. Lowering the scat- 
tering cross section of the grains, and lowering the degree of 
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Figure 5. The equatorial optical depth for scattering, absorption, 
and total extinction through the torus as a function of wave- 
length. 



Figure 6. The maximum degree of polarization of single scatter- 
ing as a function of wavelength. The maximum value is usually 
around a scattering angle of 90 degrees. The values derived from 
the observations have been found by scaling the model images 
to match the observations. The scaling factor directly gives the 
scaling that has to be applied to the model particles. 



polarization of the grains. Lowering the scattering cross sec- 
tion is difficult because it implies lowering the overall optical 
depth. This means that the torus has a lower optical depth, al- 
ready at 0.58 Jim which is inconsistent with the decreased equa- 
torial brightness seen in the ExPo images. Lowering the degree 
of polarization seems more likely. In the model grains we use, 
the degree of polarization is roughly constant over the visible 
part of the spectrum. The observations seem to indicate that the 
degree of polarization at wavelengths shorter than 0.7 jim is sig- 
nificantly lower than we predict, up to almost a factor of five at 
0.58 Jim. This can be eff'ects of surface roughness or internal ir- 
regularities of the grains. The grains we model have a grey polar- 
ization color, i.e. the degree of polarization is unchanged when 
going from 580 to 760 nm wavelength. In reality the color of po- 
larization is likely red, i.e. the degree of polarization is lower at 
shorter wavelengths (like for example is the case for cometary 
dust Kolokolova & Kimura|[2010| ). In Fig. [6] we show the max- 
imum degree of polarization, usually occurring around 90 de- 
grees scattering, as a function of wavelength. We show both the 
values for the model particles as well as the values inferred from 
the mismatch between model and ExPo images. These observa- 
tionally derived intrinsic degrees of polarization of the dust par- 
ticles are derived by equating the model images and the obser- 
vations. These values have large uncertainties, but seem to be in 
line with the steady decline with decreasing wavelength as seen 
wavelengths longwards of L6 jim. However, the model curve in- 
creases with decreasing wavelength shortwards of 1.6yL/m. This 
is likely due to an increased contribution from surface scatter- 
ing at short wavelengths. Surface scattering off' a smooth surface 
causes high degrees of polarization at the Brewster angle. This 
is not observed for scattering off' rough surfaces. We therefore 
tentatively attribute the diff'erence between model and observed 
degree of polarization in the optical to surface roughness of the 
grains. Interestingly, at wavelengths longwards of L6 jim the de- 
gree of polarization predicted by our model is in line with the ob- 
servations (see Section 4.4 ). This is what one would expect from 
this argumentation, since at long wavelengths surface eff'ects are 
negligable. Detailed modeling of this is computationally chal- 
lenging and beyond the scope of this paper. 



4.3. HST image 

|Ueta et aL] ( |2000| ) present an image taken with the Hubble Space 
Telescope (HST) of HD 161796. This image displays a strong 
centrally peaked object with no clearly resolved outer structure. 
For an optically thin shell, one would expect to see the limb 
brightened inner edge of the shell clearly resolved. However, 
since the shell is optically thick at visible wavelengths, this ef- 
fects is not observed. The polarized ffux images enhance the 
outer edges of the shell in the image, since the forward scattered 
light is unpolarized, and thus filtered out. The projected edges 
of the shell scatter at approximately 90 degrees scattering angle, 
and are thus highly polarized. This makes detecting the edge of 
the shell possible in polarized ffux (our ExPo images), but not 
show up in scattered light images of total intensity (HST im- 
age). Indeed, a simulated HST observation of our mo del shows 
a morp hology consistent with the image presented in ,Ueta et al. 
pOOO) (see also the middle panel of Fig.[TO| which is discussed 
in section [5]). 

4.4. Near infrared images 



Gledhill et al"] ( [20QT 



have obtained polarimetric images of 
HD 161796 in the near infrared using the 3.8 m United Kingdom 
Infrared Telescope (UKIRT) on Mauna Kea. The obtained im- 
ages are in the J and K band and they ffnd that there are two 
bright blobs located east and west of the central star. This struc- 
ture is more pronounced in the longer wavelength K-band image 
than in the J-band image. They report degrees of polarization up 
to ~ 20 %. We constructed model images at these wavelength, 
convolved them with a psf from a 3.8 m mirror and smeared them 
with a seeing of 0.6'', similar to the resolution they claim. Also 
we added some photon noise to simulate the observations. The 
resulting im ages are shown in Fig. [7] Comparing these images 
to Fig. 3 in |Gledhill et al.| f2001) we see that we have a very 
nice match. Indeed, at these wavelengths the equatorial struc- 
ture dominates the scattered light. Also, the modeled degree of 
polarization at these wavelengths matches the observations. 
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Figure?. Model images of the polarized flux in the J (upper 
panel) and K (lower panel) bands. To be compared with obser- 
vations from |Gledhill et al.| ( |2QQl| ). The colormap is the same as 

Fig. a 



4.5. Mid infrared images 

Mid infrared images of HD 1617 96 have been obta i ned by 
1M4\ and later by [Gledhin & ^tes| ( [2003] ). 



Skinner et al. 



Around 10 //m these images display two bright blobs east and 
west of the central star, while the central star itself is not promi- 
nent in the image. A reconstructed image at 12.5//m convolved 
with a psf of 0.5'', similar to the observations, does not fully re- 
produce this behavior. We find that we still see a contribution 
from the central star (see Fig. [8|. We have put significant ef- 
fort into solving this problem, but without success. Increasing 
the contrast between the equatorial structure and the star re- 
quires additional absorption in the infrared. However, we find 
that when we do this the total reprocessed luminosity goes up, 
so the SED does not fit anymore. A way to resolve this would be 
reducing the optical absorption while keeping the infrared ex- 
tinction the same. We can do this by removing the iron from the 
silicate grains. However, we find that in this case the silicates be- 
come too cold and we have to put them closer in to fit the SED. 
However, this is in contrast with the images that actually con- 
strain the inner radius of the dust grains very well. In addition, 
increasing the optical depth at lOyum causes the bright blobs to 
be less distinctive, and the fit to the SED becomes much worse. 
We conclude that at this moment we do not have a solution for 
this apparent discrepancy. 

We n o te here that the total infrared flux reported by 
& Yates I (|2003|) from the images is significantly hig 
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Figure 8. Model image at 12.5 yum. 

those obtained from the ISO spectrum or the fluxes measured 
by |Hrivnak et al.| ( |1989) and |Meixner et al.| (|1999|). If we ad- 
just the model to match the flux derived by Gledhill & Yates 
(2003 ) we do find a solution where the bright torus dominates 
over the light from the central resolution element. However, the 
higher flux levels reported there cannot be explained by the mild 
variations of the central source. Also, we note that we did not 
succeed in reproducing the fit to the infrared images obtained by 



Gledhill & Yates ( 2003) or Meixner et al. ( 2002). Unfortunately, 



the model images in these papers are only presented in contour 
plots and not in detailed grey- scale plots, which makes a detailed 
comparison difficult. At this point we must conclude this is still 
an unsolved issue. 



5. Mass loss history 



GledhiU 



We find a slightly smaller inner radius than [Hoogzaad et al.| 
(2002), so we find that the mass loss of HD 161796 sloped 
around 285 years ago in stead of 430 years as computed by 
[Hoogzaad et al.| ( [2QQ2| ). When we assume a constant outflow ve- 
locity (we take 15 km/s), we can compute the mass loss history 
of the object from our derived density structure. This is shown 
in Fig. [21 

It is interesting to see what the derived density structure of 
the circumstellar shell would have looked like in the past. We can 
easily compute backwards, assuming a constant outflow velocity 
of the material, to take snapshots of the nebula at difl'erent times. 
Also, we can predict what the nebula will look like in a thou- 
sand years. Model images are shown in Fig. 10 for A = 0.75 yum, 
and convolved with a PSF typical for a HST type telescope. We 
show here the image from 250 years ago since the image com- 
puted for 285 years ago, right after the mass loss stopped, is 
completely obscured at visible wavelengths. It is clearly seen 
that in the evolution from 250 years ago, when the star was still 
a heavily obscured OH/IR type star, to 1000 years from now, 
the image morphology changes significantly. This is an eff'ect 
of the changing optical depth through the circumstellar enve- 
lope. In the image close after the mass loss stopped, the optical 
depth through the shell is still very high. This results in the bipo- 
lar structure observed in the leftmost panel of Fig. [TO] At the 



ler than 



present day, the optical depth through the shell is much less, and 
the central bright component dominates. In the image computed 
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Figure 9. Mass loss history of HD 161796 as derived from the 
density of the model. The total mass loss rat e is plotted using a 
gas/dust ratio of 270 poogzaad et al.||20Q2|). 



fort into solving this problem, without success. This is an issue 
which should be addressed in future detailed modeling of this 
object. 

We find that the degree of polarization observed is signif- 
icantly lower than what is computed using a simplified dust 
model. Model particles show an unrealistically high degree of 
polarization caused by scattering off' an infinitely smooth sur- 
face. Therefore, we attribute the diff'erence between observation 
and model to surface roughness and small scale irregularities of 
the dust grains surrounding HD 161796. 

Finally, we find that the central star is significantly hotter 
than often assumed. This is mainly caused by the fact that the 
star is obscured almost completely by the equatorial structure. 
This causes a change of the spectrum of the star which makes it 
appear cooler than it actually is. This is similar to the recent anal- 
ysis by Kipper ( 2007 ). For this high temperature to be consistent 
with the observed SED, it is needed to significantly redden the 
source. Kipper] ( [2007] ) suggest that this reddening is interstellar, 
requiring a very large distance, while we show that the redden- 
ing is circumstellar, so we do not need to increase the distance 
to the source. 



for 1000 years in the future (rightmost panel in Fig. 10 ) the enve- 
lope has reached a distance from the central star far enough to be 
resolved. Also, the optical depth through the shell has decreased 
such that now the arcs to the east and west dominate in stead of 
the bipolar structure seen at early times. 

6. Conclusions 

We have imaged HD 161796 in the optical using imaging po- 
larimetry in various filters. We show that the structure seen varies 
from dominated by the low density polar outffow at short wave- 
length to a structure dominated by the high density equatorial re- 
gions at long wavelengths. This is easily explained by an equato- 
rial density enhancement which is optically thick at short wave- 
lengths and changes to optically thin at longer wavelengths. We 
show that the point where this change in geometry takes place for 
the dust shell around HD 161796 lies around 1 yum, where the op- 
tical depth through the equator is around unity. Ueta et al. ( 2000] ) 
classified the two diff'erent appearances of proto-planetary nebu- 
lae as SOLE (optically thin torus) and DUPLEX (optically thick 
torus). We show that for HD 161796 this classification is wave- 
length dependent showing that indeed these are two diff'erent 
manifestations of the same physical structure. 

We present a model for the circumstellar shell which is con- 
sistent with most observational constraints. We find that in the 
optical the integrated ffux is dominated by scattered light, while 
the central star is obscured from sight by the ffattened, equato- 
rial structure. This implies that the reddening of the central star 
is predominantly circumstellar. In order to explain our polari- 
metric images we need to have dust grains present in the torus 
with a relatively high albedo of ~ 0.84. We achieve this by using 
grains with moderate sizes. When the grains are too large, the 
albedo becomes too high, when they are too small, the albedo 
becomes too low. Given the assumptions we made for the com- 
position and shape distribution of the grains we can constrain the 
size distribution accurately to range from 0.005 - 0.5 jim with a 
power law distribution given by n(a)da oc a~^-^ . 

Our model is consistent with the SED and with the images 
available from 0.58 - 2.2yum. The model is roughly consistent 
with the observed mid infrared images with the important ex- 
ception that the central star remains visible in our model, while 
it is not seen in the observed images. We have put significant ef- 
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Figure 10. Model images of the total intensity as a function of time, assuming a constant outflow velocity of 15 km/s. The model 
images are copmputed at A = 0.75 yum, and convolved with a point spread function representative for an HST type telescope. The 
left image is computed for 250 years ago, the middle image is for the present day, and the rightmost image is computed for 1000 
years from now. 
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